Graph isomorphism and adiabatic quantum computing 
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In the Graph Isomorphism problem two iV-vertex graphs G and G' are given and the task is to 
determine whether there exists a permutation of the vertices of G that preserves adjacency and 
transforms G — > G' . If yes, then G and G' are said to be isomorphic; otherwise they are non- 
isomorphic. The GI problem is an important problem in computer science and is thought to be of 
comparable difficulty to integer factorization. In this paper we present a quantum algorithm that 
solves arbitrary instances of GI and can also determine all automorphisms of a given graph. We show 
how the GI problem can be converted to a combinatorial optimization problem that can be solved 
using adiabatic quantum evolution. We numerically simulate the algorithm's quantum dynamics and 
show that it correctly: (i) distinguishes non- isomorphic graphs; (ii) recognizes isomorphic graphs; 
and (iii) finds all automorphisms of a given graph G. We then discuss the GI quantum algorithm's 
experimental implementation, and close by showing how it can be leveraged to give a quantum 
algorithm that solves arbitrary instances of the NP-Complete Sub-Graph Isomorphism problem. 

PACS numbers: 03.67.Ac,02.10.Ox,89.75.Hc 



I. INTRODUCTION 

An instance of the Graph Isomorphism (GI) problem 
is specified by two iV-vertex graphs G and G' and the 
challenge is to determine whether there exists a permu- 
tation of the vertices of G that preserves adjacency and 
transforms G — > G' . When such a permutation exists, 
the graphs are said to be isomorphic; otherwise they are 
non-isomorphic. GI has been heavily studied in com- 
puter science [l|. Polynomial classical algorithms exist 
for special cases of GI, still it has not been possible to 
prove that GI is in P. Although it is known that GI is 
in NP, it has also not been possible to prove that it is 
NP-Complete. The situation is the same for Integer Fac- 
torization (IF) — it belongs to NP, but is not known to be 
in P or to be NP-Complete. GI and IF are believed to 
be of comparable computational difficulty [2]. 

IF and GI have also been examined from the perspec- 
tive of quantum algorithms and both have been con- 
nected to the hidden subgroup problem (HSP) [3]. For 
IF the hidden subgroup is contained in an abelian parent 
group (Z* — group of units modulo n) , while for GI the 
parent group is non-abelian (S n — symmetric group on 
n elements). While Fourier sampling allows the abelian 
HSP to be solved efficiently [4], strong Fourier sampling 
does not allow an efficient solution of the non-abelian 
HSP over S n [a|. At this time an efficient quantum algo- 
rithm for GI is not known. 

A number of researchers have considered using the 
dynamics of physical systems to solve instances of GI. 
Starting from physically motivated conjectures, these ap- 
proaches embed the structure of the graphs appearing in 
the GI instance into the Hamiltonian that drives the sys- 
tem dynamics. In Ref. 6| the systems considered were 
classical, while Refs. B-tllj worked with quantum sys- 
tems. See Refs. [10|, LLll for a summary of the current 
status of this approach to GI. 



In this paper we present a quantum algorithm that 
solves arbitrary instances of GI. The algorithm can also 
be used to determine all automorphisms of a given graph. 
The GI quantum algorithm is constructed by first con- 
verting an instance of GI into an instance of a combi- 
natorial optimization problem whose cost function has a 
zero minimum value when the pair of graphs in the GI 
instance are isomorphic, and is positive when the pair are 
non-isomorphic. The construction of the GI quantum al- 
gorithm is completed by showing how the combinatorial 
optimization problem can be solved using adiabatic quan- 
tum evolution. To test the efficiacy of this GI quantum 
algorithm we numerically simulated its Schrodinger dy- 
namics. The simulation results show that it can correctly: 
(i) distinguish pairs of non-isomorphic graphs; (ii) recog- 
nize pairs of isomorphic graphs; and (iii) find all automor- 
phisms of a given graph. We also discuss the experimen- 
tal implementation of the GI algorithm, and show how it 
can be leveraged to give a quantum algorithm that solves 
arbitrary instances of the (NP-Complete) SubGraph Iso- 
morphism problem. 

The structure of this paper is as follows. In Section [TT1 
we give a careful presentation of the GI problem, and 
show how an instance of GI can be converted to an in- 
stance of a combinatorial optimization problem whose 
solution (Section lIII[) can be found using adiabatic quan- 
tum evolution. To test the performance of the GI quan- 
tum algorithm introduced in Section [ill! we numerically 
simulated its Schrodinger dynamics and the results of 
that simulation are presented in Section [TV] In Section IVl 
we describe the experimental implementation of the GI 
algorithm, and in Section IVll we show how it can be used 
to give a quantum algorithm that solves arbitrary in- 
stances of the NP-Complete problem known as SubGraph 
Isomorphism. Finally, we summarize our results in Sec- 
tion En 



II. GRAPH ISOMORPHISM PROBLEM 

In this Section we introduce the Graph Isomorphism 
(GI) problem and show how an instance of GI can be con- 
verted into an instance of a combinatorial optimization 
problem (COP) whose cost function has zero (non-zero) 
minimum value when the pair of graphs being studied 
are isomorphic (non- isomorphic) . 



A. Graphs and graph isomorphism 

A graph G is specified by a set of vertices V and a 
set of edges E. We focus on simple graphs in which an 
edge only connects distinct vertices, and the edges are 
undirected. The order of G is defined to be the number 
of vertices contained in V, and two vertices are said to be 
adjacent if they are connected by an edge. If x and y are 
adjacent, we say that y is a neighbor of x, and vice versa. 
The degree d(x) of a vertex x is equal to the number of 
vertices that are adjacent to x. The degree sequence of a 
graph lists the degree of each vertex in the graph ordered 
from largest degree to smallest. A graph G of order N 
can also be specified by its adjacency matrix A which is 
an N x N matrix whose matrix clement 



'■'•] 



1(0) if 

the vertices i and j are (are not) adjacent. For simple 
graphs a,i t i = and ai.j = aj.i since edges only connect 
distinct vertices and are undirected. 

Two graphs G and G' are said to be isomorphic if there 
is a one-to-one correspondence 7r between the vertex sets 
V and V such that two vertices x and y are adjacent in 
G if and only if their images n x and n y are adjacent in 
G' . The graphs G and G' are non-isomorphic if no such 
7r exists. Since no one-to-one correspondence it can exist 
when the number of vertices in G and G' are different, 
graphs with unequal orders are always non-isomorphic. It 
can also be shown [12J | that if two graphs are isomorphic, 
they must have identical degree sequences. 

We can also describe graph isomorphism in terms of 
the adjacency matrices A and A' of the graphs G and 
G', respectively. The graphs are isomorphic if and only 
if there exists a permutation matrix a of the vertices of 
G that satisfies 



A' = aAa T , 



(1) 



where o~ T is the transpose of a. It is straight-forward to 
show that if G and G" are isomorphic, then a permuta- 
tion matrix a exists that satisfies Eq. §T§. To prove the 
"only if" statement, note that the i-j matrix element of 
the RHS is A ff- . a .. Eq. (fTJ) is thus satisfied when a per- 
mutation matrix a exists such that A', „■ = A„. „.. Thus 
Eq. (JTJ) is simply the condition that a preserve adjacency. 
Since a permutation is a one-to-one correpondence, the 
existence of a permutation matrix a satisfying Eq. ([1]) 
implies G and G" are isomorphic. Thus, the existence of 
a permutation matrix a satisfying Eq. |T]) is an equivalent 
way to define graph isomorphism. 



The Graph Isomorphism (GI) problem is to determine 
whether two given graphs G and G" are isomorphic. The 
problem is only non-trivial when G and G' have the same 
order and so we focus on that case in this paper. 



B. Permutations, binary strings, and linear maps 

A permutation it of a finite set S = {0, . . . , N — 1} is 
a one-to-one correspondence from 5 — ?• <S which sends 
i — > iTi such that -Ki £ S, and -Ki =/= ttj for i =/= j- The 
permutation 7r can be written 







I 



N -1 

7TJV-1 



(2) 



where column i indicates that w sends i — > 7Tj. Since 
the top row on the RHS of Eq. ^ is the same for all 
permutations, all the information about 7r is contained 
in the bottom row. Thus we can map a permutation -k 
into an integer string P(ir) = ttq • • -ttn-i, with iri £ S 
and iTi ^ ttj for i ^ j . 

For reasons that will become clear in Section III CI we 
want to convert the integer string P(n) = ttq ■ ■ ■ t^n-i 
into a binary string p^. This can be done by replacing 
each -Ki in P(tt) by the unique binary string formed from 
the coefficients appearing in its binary decomposition 



U-i 

n i = 5Z ni >3 ( 2 ^ 
3=0 



(3) 



Here U = [log 2 N~\ . Thus the integer string P(ir) is trans- 
formed to the binary string 

Ptt = (tTo.O ' ' ' 7To,U-l) • • ■ (7TJV-1,0 ' ' ' 1fN-l,U-l) i (4) 

where ir.ij £ {0, 1}. The binary string p n has length NU, 
where 



N <2 U = M+1. 



(5) 



Thus we can identify a permutation ir with the binary 
string pn in Eq. (g|). 

Let T-L be the Hamming space of binary strings of length 
NU. This space contains 2 NU strings, and we have just 
seen that AH of these strings p^ encode permutations 
7T. Our last task is to define a mapping from H to the 
space oi N x N matrices a with binary matrix elements 
o~i,j — 0, 1. The mapping is constructed as follows: 

1. Let Sb = so • ■ ■ snu^i be a binary string in T-L. We 
parse s& into TV substrings of length U as follows: 

Sb = {so ■ ■ ■ SU-l) (SU ■ • • S 2 U-l) ■ ■ ■ (s(Ar_i)(7 • ■ • SOT-l) • 

(6) 

2. For each substring sm ■ ■ ■ s^+iw-i! construct the 
integer 



(7-1 

E 

j=0 



Siu+j (2) J < 2 U - 1 = M. 



(7) 



Finally, introduce the integer string s = 
So ■ ■ ■ sjv-i, and define the N x N matrix a(s) to 
have matrix elements 



,•(*) 



0, iisj>N-l 
&„, if 0<*,- < JV-1, 



(8) 



where i, j 6 <S, and 6 x ,y is the Kronecker delta. 

Note that when the binary string si, corresponds to a per- 
mutation, the matrix cr(s) is a permutation matrix since 
the Si formed in step 2 will obey < Sj < N ~ 1 and 
Si y^= Sj for i ^ j. In this case, if A is the adjacency 
matrix for a graph G, then A' — a(s)Acr T (s) will be the 
adjacency matrix for a graph G" isomorphic to G. On 
the other hand, if s& does not correspond to a permuta- 
tion, then the adjacency matrix A' = ct(s)j4<t t (s) must 
correspond to a graph G" which is not isomorphic to G. 
The result of our development so far is the establish- 
ment of a map from binary strings of length NU to N x N 
matrices (viz. linear maps) with binary matrix elements. 
When the string is (is not) a permutation, the matrix 
produced is (is not) a permutation matrix. 



C. Graph isomorphism and combinatorial 
optimization 



Here ||M||, is the L^-norm of M. In the numerical sim- 
ulations discussed in Section IIV[ the ii-norm is used, 
though any Li-norm would be acceptable. Thus, when 
G and G' are isomorphic, C*3(s) = 0, and a(s) is the per- 
mutation of vertices of G that maps G — > G 1 '. Putting all 
these remarks together, we see that if C(s) — for some 
integer string s, then G and G' are isomorphic and <r(s) 
is the permutation that connects them. On the other 
hand, if C(s) > for all strings s, then G and G' are 
non-isomorphic. 

We have thus converted an instance of GI into an in- 
stance of the following COP: 

Graph Isomorphism COP: Given the TV-vertex 
graphs G and G' and the associated cost function 
C(s) defined above, find an integer string s* that 
minimizes C(s). 

By construction: (i) C(s*) = if and only if G and G' are 
isomorphic and <r(s*) is the permutation matrix mapping 
G ->■ G"; and (ii) C(s Sf ) > if and only if G and G" are 
non-isomorphic. 

Before moving on, notice that if G = G', then C(s*) = 
since G is certainly isomorphic to itself. In this case 
cr(s*) is an automorphism of G. We shall see that the GI 
quantum algorithm to be introduced in Section IIIII can 
thus be used to find the automorphism group of a graph. 



As seen above, an instance of GI is specified by a pair of 
graphs G and G" (or equivalently, by a pair of adjacency 
matrices A and A'). Here we show how a GI instance 
can be transformed into an instance of a combinatorial 
optimization problem (COP) whose cost function has a 
minimum value of zero if and only if G and G' are iso- 
morphic. 

The search space for the COP is the Hamming space 
"H of binary strings s& of length NU which are associated 
with the integer strings s and matrices <j(s) introduced in 
Section Hi Bl The COP cost function C(s) contains three 
contributions 



C(s) = C 1 (s)+C 2 (s)+C 3 (s). 



(9) 



The first two terms on the RHS penalize integer strings 
s = sq' • • Sff—i whose associated matrix a(s) is not a 
permutation matrix, 



C 1 (s) 



C 2 (s) 



N-l M 

E E *«,« 

i=0 a=N 
N-2 N-l 

E E *«- 

i=0 j=i+l 



(10) 



(11) 



where S XlV is the Kronecker delta. We see that Ci(s) > 
when Si > N — 1 for some i, and C2{s) > when Sj = Sj 
for some i jt j. Thus Ci(s) + C?2(s) = if and only if 
cr(s) is a permutation matrix. The third term Csls) adds 
a penalty when a(s)Aa T (s) ^ A': 



III. ADIABATIC QUANTUM ALGORITHM 
FOR GRAPH ISOMORPHISM 



The adiabatic quantum evolution (AQE) algorithm 
13] exploits the adiabatic dynamics of a quantum system 
to solve COP. The AQE algorithm uses the optimization 
problem cost function to define a problem Hamiltonian 
Hp whose ground-state subspace encodes all problem so- 
lutions. The algorithm evolves the state of an L-qubit 
register from the ground-state of an initial Hamiltonian 
Hi to the ground-state of Hp with probability approach- 
ing 1 in the adiabatic limit. An appropriate measurement 
at the end of the adiabatic evolution yields a solution of 
the optimization problem almost certainly. The time- 
dependent Hamiltonian H(t) for global AQE is 



H(t) 



1 



H, 



Hp, 



(13) 



where T is the algorithm runtime, and adiabatic dynam- 
ics corresponds to T — > oo. 

To map the GI COP onto an adiabatic quantum com- 
putation, we begin by promoting the binary strings s& to 
computational basis states (CBS) |s&). Thus each bit in 
Sb is promoted to a qubit so that the quantum register 
contains L = NU = N\\og 2 N~\ qubits. The CBS are 
defined to be the 2 L eigenstates of cr° <8> ■ • • ® erf' -1 . The 
problem Hamiltonian Hp is defined to be diagonal in the 
CBS with eigenvalue G(s), where s is the integer string 
associated with Sb- 



C 3 (s)=\\a(s)Aa T (s)-A'\\ l 



(12) 



H P \s b ) = C(s)\s b ). 



(14) 



Note that (see Section III C|) the ground-state energy of 
Hp will be zero if and only if the graphs G and G' are 
isomorphic. We will discuss the experimental realization 
of Hp in Section [V] The initial Hamiltonian Hi is chosen 
to be 



H, 



L-l 

1=0 



oi) 



(15) 



where I and a l x are the identity and x-Pauli operator for 
qubit I, respectively. The ground-state of Hi is the easily 
constructed uniform superposition of CBS. 

The quantum algorithm for GI begins by preparing the 
L qubit register in the ground-state of Hi and then driv- 
ing the qubit register dynamics using the time-dependent 
Hamiltonian H{t). At the end of the evolution the qubits 
are measured in the computational basis. The outcome 
is the bit string si so that the final state of the register is 
\sl) and its energy is C(s*), where s* is the integer string 
derived from s£. In the adiabatic limit, C(s*) will be the 
ground-state energy, and if C(s*) = ( > 0) the algo- 
rithm decides G and G' are isomorphic (non-isomorphic) . 
Note that any real application of AQE will only be ap- 
proximately adiabatic. Thus the probability that the fi- 
nal energy C(s*) will be the ground-state energy will be 
1 — e. In this case the GI quantum algorithm must be 
run k ~ C(ln(l — 6)/ lne) times so that, with probability 
5 > 1 — e, at least one of the measurements will return 
the ground-state energy. We can make S arbitrarily close 
to 1 by choosing k sufficiently large. 



IV. NUMERICAL SIMULATION OF 
ADIABATIC QUANTUM ALGORITHM 

In this Section we numerically simulate the dynamics 
of the GI adiabatic quantum algorithm ( AQA) . Because 
the GI AQA uses N\log 2 N~\ qubits, and these simula- 
tions were carried out using a classical digital computer, 
we were limited to GI instances involving graphs with 
order N < 7. Although we would like to have examined 
larger graphs, this simply was not practical. Note that 
the N = 7 simulations use 21 qubits. These simulations 
are at the upper limit of 20-22 qubits at which simula- 
tion of the full adiabatic Schrodinger dynamics is feasible 
[lj-[l6| • To simulate a GI instance with graphs of order 
N = 8 would require a 24 qubit simulation which is well 
beyond what can be done practically. The protocol for 
the simulations presented here follows Refs. |L4H17| . 

In Section llV Al we present simulation results for sim- 
ple examples of isomorphic and non-isomorphic graphs. 
These examples allow us to illustrate the analysis of the 
simulation results in a simple setting. Section IIVBI then 
presents our simulation results for non-isomorphic in- 
stances of: (i) iso-spectral graphs; and (ii) strongly regu- 
lar graphs. Finally, Section TlV CI considers GI instances 
where G' — G. Clearly, all such instances correspond 
to isomorphic graphs since the identity permutation will 



always map G —> G and preserve adjacency. The situ- 
ation is more interesting when G has symmetries which 
allow non-trivial permutations as graph isomorphisms. 
These self-isomorphisms are referred to as graph auto- 
morphisms, and they form a group known as the auto- 
morphism group Aut(G) of G. In this final subsection we 
use the GI AQA to find Aut(G) for a number of graphs. 
In all GI instances considered in this Section, the GI 
AQA correctly: (i) distinguished non-isomorphic pairs 
of graphs; (ii) recognized isomorphic pairs of graphs; 
and (iii) determined the automorphism group of a given 
graph. 



A. Illustrative examples 

Here we present the results of a numerical simulation 
of the GI AQA applied to two simple GI instances. In 
Section IIV All (Section IIV A2[) we examine an instance 
of two non-isomorphic (isomorphic) graphs. For the iso- 
morphic instance we also present the permutations found 
by the GI AQA that transforms G into G' while preserv- 
ing adjacency. 



1. Non-isomorphic graphs 

Here we use the GI AQA to examine a GI instance in 
which the two graphs G and G" are non-isomorphic. The 
two graphs are shown in Figure [TJ Each graph contains 




G 



FIG. 1: Two non-isomorphic 4-vertex graphs G and G' . 

4 vertices and 5 edges, however they are non-isomorphic. 
Examining Figure Q] we see that the degree sequence for 
G is {2,2,2,2}, while that for G' is {3,2,2,1}. Since 
these degree sequences are different we know that G and 
G' are non-isomorphic. Finally, the adjacency matrices 
A and A' for G and G", respectively, are: 



.4 



/0 1 1 
10 10 
10 1 

\1 1 



A 1 




(16) 



The GI AQA finds a non-zero final ground-state energy 
E gs = 4, and so correctly identifies these two graphs 
as non-isomorphic. It also finds that the final ground- 
state subspace has a degeneracy of 16. The linear maps 
associated with the 16 CBS that span this subspace give 
rise to the lowest cost linear maps of G relative to G' . As 
these graphs are not isomorphic, these lowest cost maps 
have little inherent interest and so we do not list them. 



2. Isomorphic graphs 

Here we examine the case of two isomorphic graphs G 
and G' which are shown in Figure [5J Each graph con- 





FIG. 2: Two isomorphic 4-vertex graphs G and G' . 

tains 4 vertices and 5 edges, and both graphs have degree 
sequence {3, 3, 2, 2}. By inspection of Figure^ the asso- 
ciated adjacency matrices are, respectively, 



A = 




A' = 




(17) 



For these two graphs, the GI AQA finds a vanishing final 
ground-state energy E gs = and so recognizes G and G" 
as isomorphic. It also finds that the final ground-state 
subspace has a degeneracy of 4. The four CBS that span 
this subspace give four binary strings Sb- These four bi- 
nary strings in turn generate four integer strings s which 
are, respectively, the bottom row of four permutations 
(see Eq. (|2J)). Thus, not only does the GI AQA recognize 
G and G' as isomorphic, but it also returns the 4 graph 
isomorphisms ~k\ , . . . , n^ that transform G — > G' while 
preserving adjacency: 



71"! 



7T3 



12 3 

2 3 1 

12 3 

3 10 2 



7T2 = 



7T 4 



12 3 

3 2 1 

12 3 

13 2 



(18) 



We will next explicitly show that tti is a graph isomor- 
phism; the reader can easily check that the remaining 3 
permutations are also graph isomorphisms. 

The permutation matrix a\ associated with ~k\ is 



<7l = 



/ 1 

1 

10 

\0 1 



(19) 



Under tti, G is transformed to tt\(G) which is shown 
in Figure [3) It is clear from Figure [3] that vertices x 
and y are adjacent in G if and only if iri tX and tti^ v are 
adjacent in 7Ti(G). The adjacency matrix for 7Ti(G) is 
<j\Aa\ which is easily shown to be 



oxAo[ = 




(20) 





FIG. 3: Transformation of G produced by the permutation 

7Tl. 



This agrees with the adjacency shown in 7Ti(G) in Fig- 
ure|3] Comparison with Eq. (fT7|) shows that o\Ao\ = A'. 
Thus the permutation -k\ does map G into G' and pre- 
serve adjacency and so establishes that G and G' are 
isomorphic. 

Finally, note that G and G' are connected by exactly 
4 graph isomorphisms. Examination of Figure [5] shows 
that the degree-3 vertices in G are vertices and 2, while 
in G" they are vertices and 3. Since the degree of a 
vertex is preserved by a graph isomorphism [12j, vertex 

in G must be mapped to vertex or 3 in G'. Then 
vertex 2 (in G) must be mapped to vertex 3 or (in G'), 
respectively. This then forces vertex 1 to map to vertex 

1 or 2, and vertex 3 to map to vertex 2 or 1, respectively. 
Thus only 4 graph isomorphism are possible and these 
are exactly the 4 graph isomorphisms found by the GI 
AQA which appear in Eq. (Ql 



B. Non-isomorphic graphs 

In this subsection we present GI instances involving 
pairs of non-isomorphic graphs. In Section [IV B 1 1 we ex- 
amine two instances of isospectral graphs; and in Sec- 
tion [TvE2] we look at three instances of strongly regular 
graphs. We shall see that the GI AQA correctly distin- 
guishes all graph pairs as non-isomorphic. 



1. Iso-spectral graphs 

The spectrum of a graph is the set containing all 
the eigenvalues of its adjacency matrix. Two graphs 
are iso-spectral if they have identical spectra. Non- 
isomorphic iso-spectral graphs are believed to be difficult 
to distinguish [8|, Q • Here we test the GI AQA on pairs of 
non-isomorphic iso-spectral graphs. The non-isomorphic 
pairs of iso-spectral graphs examined here appear in 
Ref. 0. 

N = 5 : It is known that no pair of graphs with less 
than 5 vertices is iso-spectral [18J . Figure 2] shows a pair 
of graphs G and G' with 5 vertices which are isospec- 
tral and yet are non-isomorphic. Although both have 
5 vertices and 4 edges, the degree sequence of G is 
{2, 2, 2, 2, 0}, while that of G' is {4, 1, 1, 1, 1}. Since they 
have different degree sequences, it follows that G and G' 





FIG. 4: Two non-isomorphic 5-vertex iso-spectral graphs G 
and G'. 



are non-isomorphic. The adjacency matrices A and A' 
for G and G", respectively, are: 



.4 



/0 1 1 0\ 

10100 1 
10 10 
10 10 

\o o o o/ 



A' 



/0 1 1 1 1 \ 

10 0* 

10 

10 

V i o o o o/ 



(21) 



It is straight-forward to show that A and A' have the 
same characteristic polynomial P{\) = A 5 — 4A 3 and 
so G and G' are iso-spectral. Numerical simulation 
of the dynamics of the GI AQA finds a non-zero final 
ground-state energy E gs = 5, and so the GI AQA 
correctly distinguishes G and G' as non-isomorphic 
graphs. 

N = 6 : The pair of graphs G and G' in Figure [5] arc 

0. A 2 

4 5 





FIG. 5: Two non-isomorphic 6-vertex iso-spectral graphs G 
and G". 

iso-spectral and non-isomorphic. To see this, note that 
although both graphs have 6 vertices and 7 edges, the 
degree sequence of G is {5, 2, 2, 2, 2,1}, while that of G' is 
{3, 3, 3, 3, 1, 1}. Since their degree sequences are different, 
G and G" are non-isomorphic. The adjacency matrices A 
and A 1 of G and G', respectively, are: 



.4 



/0 1 1\ 


' 1 1 


1 


11 


10 1 


V i i i i i o/ 



A' 



/0 1 0\ 


10 110 


10 110 


110 10 


110 1 


\0 1 0/ 



(22) 
Both A and A' have the same characteristic polynomial 
P(\) = A 6 - 7A 4 - 4A 3 + 7A 2 + 4A - 1 which establishes 



that G and G' are iso-spectral. Numerical simulation 
of the dynamics of the GI AQA finds a non-zero final 
ground-state energy E gs — 7, and so the GI AQA cor- 
rectly distinguishes G and G' as non-isomorphic graphs. 



2. Strongly regular graphs 



As we saw in Section Hi A[ the degree d(x) of a vertex 
x is equal to the number of vertices that are adjacent to 
x. A graph G is said to be k-regular if, for all vertices x, 
the degree d(x) = k. A graph is said to be regular if it is 
/c-regular for some value of k. Finally, a strongly regular 
graph is a graph with v vertices that is fc-regular, and 
for which: (i) any two adjacent vertices have A common 
neighbors; and (ii) any two non-adjacent vertices have fj, 
common neighbors. In this subsection we apply the GI 
AQA to pairs of non-isomorphic strongly regular graphs. 

N = 4 : In Figure [6] we show two strongly regular 4- 






FIG. 6: Two 4-vertex non-isomorphic strongly regular graphs 
G and G". 

vertex graphs G and G'. The parameters for G are v = 4, 
k = 2, A = 0, and fi = 2; while for G' they are v = 4, 
k = 3, A = 2, and /i = 0. It is clear that G and G' are 
non-isomorphic since they contain an unequal number 
of edges and different degree sequences. The adjacency 
matrices A and A' for G and G' are, respectively, 



A = 




A' = 




(23) 



Numerical simulation of the GI AQA dynamics finds a 
non-zero final ground-state energy E gs — 4 and so the GI 
AQA correctly distinguishes G and G' as non-isomorphic. 

N = 5 : In Figure [7] we show two strongly regular 5- 


4_/ V 1 4 // \\_ 1 





G 



FIG. 7: Two 5-vertex non-isomorphic strongly regular graphs 
G and G'. 



vertex graphs G and G' . The parameters for G are v = 5, 
fc = 2, A = 0, and fj, = 1; while for G' they are r/ = 5, 
fc = 4, A = 3, and /i = 0. It is clear that G and G' are 
non-isomorphic since they contain an unequal number 
of edges and different degree sequences. The adjacency 
matrices A and A' for G and G' are, respectively, 



A 



/ 1 1 \ 

' 1 1 ' 

10 10 

10 1 

V i o o i 0/ 



A' 



/0 1 1 1 1 \ 

' 1 1 1 1 ' 

110 11 

1110 1 

\l 1 1 1 0/ 



(24) 



Numerical simulation of the GI AQA dynamics finds 



a non-zero final ground-state energy E g 



10 and 



so the GI AQA correctly distinguishes G and G' as 
non-isomorphic. 

N = 6 : In Figure [8] we show two strongly regular 6- 

0. 





G 



G' 



FIG. 8: Two 6-vertex non-isomorphic strongly regular graphs 
G and G' . 

vertex graphs G and G'. The parameters for G are v = 6, 
k = 3, A = 0, and /i = 3; while for G' they are r/ = 6, 
fc = 4, A = 2, and n = 4. It is clear that G and G' arc 
non-isomorphic since they contain an unequal number 
of edges and different degree sequences. The adjacency 
matrices A and A 1 for G and G' are, respectively, 



.4 



Numerical simulation of the GI AQA dynamics finds 
a non-zero final ground-state energy E gs — 10 and so 
the GI AQA correctly distinguishes G and G' as non- 
isomorphic. 



C. Graph automorphisms 



/0 1 1 1 \ 
10 10 10 




/0 1 1 1 1 0\ 

401011 


10 10 1 
10 10 10 


; A' = 


110 10 1 
10 10 11 


10 10 1 




110 10 1 


V 1 1 1 0/ 




\0 1 1 1 1 0/ 



As noted in the Section IIVI introduction, we can find 
the automorphism group Aut(G) of a graph G using the 
GI AQA by considering a GI instance with G" = G. Here 
the self-isomorphisms are permutations of the vertices of 
G that map G ^ G while preserving adjacency. Since 



G is always isomorphic to itself, the final ground-state 
energy will vanish: E gs = 0. The set of CBS |sj,} that 
span the final ground-state subspace give rise to a set 
of binary strings sj, that determine the integer strings 
s = So ■ ■ ■ Sjv-i (see Section Til Bp that then determine 
the permutations 7r(s), 



tt(s) 





•so 



i 

Si 



N - 1 

SJV-I 



(26) 



which are all the elements of Aut(G). By construction, 
the order of Aut(G) is equal to the degeneracy of the 
final ground-state subspace. In this subsection we apply 
the GI AQA to the: (i) cycle graphs C4, . . . , G7; (ii) grid 
graph G2.3; and (iii) wheel graph Wr, and show that it 
correctly determines the automorphism group for all of 
these graphs. 



1. Cycle graphs 

A walk W in a, graph is an alternating sequence of 
vertices and edges xq, e±, x\, e%, ■ • ■ , e/, xi, where the 
edge &i connects Xi—\ and Xi for < i < I. A walk W 
is denoted by the sequence of vertices it traverses W = 
xo^i • • • Xi- Finally, a walk W = x$x\ ■ ■ ■ xi is a cycle if 
I > S; xq — xf. and the vertices Xi with < i < I are 
distinct from each other and xq. A cycle with n vertices 
is denoted C n , 

The automorphism group of the cycle graph C„ is the 
dihedral group D n [l2| . The order of D n is 2n, and it is 
generated by the two elements a and /3 that satisfy the 
following relations, 

a n = e; (3 2 = e; a/3 = /?a"-\ (27) 

where e is the identity element. Because a and /3 are 
generators of D n , each element g of D n can be written 
as a product of appropriate powers of a and /3: 



5 = « l /3 J 



(0 < i < n - 1 ; < j < 1). (28) 



We now use the GI AQA to find the automorphism 
group of the cycle graphs C n for 4 < n < 7. We will see 
that the GI AQA correctly determines Aut(C n ) = D n 
for these graphs. We will work out G4 in detail, and then 
give more abbreviated presentations for the remaining 
cycle graphs as their analysis is identical. 

N = 4 : The cycle graph G4 appears in Figure |9] It 
contains 4 vertices and 4 edges; has degree sequence 
{2, 2, 2, 2}; and adjacency matrix 



.4 



/0101 
10 10 
10 1 

Vi 1 



(29) 



Numerical simulation of the GI AQA applied to the 
GI instance with G = G4 and G' = G found a vanishing 



and the associated permutation matrix er* = er(3210) is 



Ca 



FIG. 9: Cycle graph C 4 . 



final ground-state energy E gs = 0. The GI AQA thus 
correctly identifies C4 as being isomorphic to itself. The 
simulation also found that the final ground-state sub- 
space is 8-fold degenerate. TableQ]lists the integer strings 



TABLE I: Automorphism group Aut{Ci) of the cycle graph 
d as found by the GI AQA. The first row lists the integer 
strings s = S0S1S2S3 determined by the labels of the eight 
CBS \sb) that span the final ground-state subspace (see text). 
Each string a determines a graph automorphism 7r(s) via 
Eq. (|26|l . The second row associates each integer string s in 
the first row with a graph automorphism n(s). It identifies 
the two strings that give rise to the graph automorphisms 
a and /3 that generate Aut{C4), and writes each graph 
automorphism n(s) as a product of an appropriate power of 
a and /3. Note that e is the identity automorphism, and the 
product notation assumes the rightmost factor acts first. 



S0S1S2S3 



tt(s) 



3012 2301 1230 0123 0321 3210 2103 1032 



P aP a 2 /3 a 6 p 



s = S0S1S2S3 that result from the eight CBS | sj,) that 
span the final ground-state subspace (see Sections III 131 
and IIIip . Each integer string s determines the bottom 
row of a permutation tt(s) (see eq. (j26|)). We explicitly 
show that 7r(s), for the integer string s = 3210, is an au- 
tomorphism of C4. The reader can repeat this analysis to 
show that the seven remaining integer strings in Table U 
also give rise to automorphisms of C4. Note that the CBS 
that span the final ground-state subspace determine all 
the graph automorphisms of C4. This follows from the 
manner in which the GI AQA is constructed since each 
automorphism of C4 must give rise to a CBS with van- 
ishing energy, and each CBS in the final ground-state 
subspace gives rise to an automorphism of C4 . These au- 
tomorphisms form a group AutiCi) and the GI AQA has 
found that the order of Aut(Ci) is 8 which is the same 
as the order of the dihedral group D4. 
The permutation 7T* = ir(s — 3210) is 



tt(3210) = 




(30) 



ct(3210) 



( {) l\ 

10 

10 

\ 1 0/ 



(31) 



Under 7T*, C4 is transformed to 71% (C4) which is shown 
in Figure [TDJ It is clear from Figure [TO] that x and y are 




FIG. 10: Transformation of C4 produced by the permutation 

it,. . 

adjacent in C4 if and only if ti*. x and ~k* iV are adjacent 
in ^(C^). Thus 7r* is a permutation of the vertices of C4 
that preserves adjacency and so is a graph automorphism 
of C4. We can also show this by demonstrating that 
o*Ao^ = A. Using Eqs. (|29|) and (|3lT) it is easy to show 
that 



tAaT 



/oioi\ 

10 10 

10 1 

\1 1 0/ 



(32) 



This agrees with the adjacency of edges in 7r*(C4) ap- 
pearing in FigureHni an d comparison with Eq. ([2"5|) shows 
that a^AaJ = A, confirming that 7r* is an automorphism 
ofC* 4 . 

We now show that Aut(C4) is isomorphic to the di- 
hedral group -D4 by showing that the automorphisms 
7r(3012) = a and 7r(0321) = f3 generate Aut(C4), and 
satisfy the generator relations (Eq. (|27[) ) for D4. The 
second row of Table |T] establishes that a and /3 are the 
generators of Aut(C4) as it shows that each element of 
Aut(C4) is a product of an appropriate power of a and 
/3, and all possible products of powers of a and (3 appear 
in that row. Now notice that 



Q 




(33) 



corresponds to a clockwise rotation of C4 by 90° (see 
Figure 111)) Thus four applications of a corresponds to 
a 360° rotation of C4 which leaves it invariant. Thus 
a 4 = e. This can also be checked by composing a with 



itself four times using Eq. (|3"31 . This establishes the first 
of the generator relations in Eq. (f2"Tf . Similarly, 



P- 




(34) 




FIG. 11: Transformation of C4 produced by the automor- 
phism a. 









/?(C 4 



FIG. 13: Cycle graph C 5 



FIG. 12: Transformation of C4 produced by the automor- 
phism p. 



corresponds to reflection of C4 about the diagonal pass- 
ing through vertices and 2 (see Figure [12]) . Thus two 
applications of j5 leaves C4 invariant, and so /3 2 = e. 
This establishes the second of the generator relations in 
Eq. (|2T|) . Finally, to show the third generator relation 
af3 — [3a 3 , we simply evaluate both sides of this relation 
and compare results. Using Table UJ we find that 



a/3 = 



(3a 3 




(35) 



(36) 



It is clear that a/3 does equal /3a 3 . Thus we have shown 
that a and (3: (i) generate Aut(Ci)\ and (ii) satisfy the 
generator relations (Eq. (|27|) for the dihedral group D4, 
and so generate a group isomorphic to D4. In summary, 
we have shown that the GI AQA found all eight graph 
automorphisms of C4, and that the group formed from 
these automorphisms is isomorphic to the dihedral group 
D4 which is the correct automorphism group for C4. 

N = 5 : The cycle graph C5 appears in Figure [131 It 
has 5 vertices and 5 edges; degree sequence {2,2,2,2,2}; 
and adjacency matrix 

/o 1 1 \ 
10 10 

A = 10 10 • (37) 

10 1 

V 1 1 0/ 



Numerical simulation of the GI AQA applied to the 
GI instance with G = C5 and G' = G found a vanish- 
ing final ground-state energy E gs = 0. The GI AQA 
thus correctly identifies C5 as being isomorphic to it- 
self. The simulation also found that the final ground- 
state subspace is 10-fold degenerate. Table [TT] lists the 



TABLE II: Automorphism group Aut{C§) of the cycle graph 
C5 as found by the GI AQA. The odd rows list the integer 
strings s = so ■ • • S4 determined by the labels of the ten CBS 
\sb) that span the final ground-state subspace (see text). 
Each string s determines a graph automorphism n(s) via 
Eq. (|26l) . Each even row associates each integer string s 
in the odd row preceding it with a graph automorphism 
7r(s). It identifies the two strings that give rise to the graph 
automorphisms a and /3 that generate Aut(C$), and writes 
each graph automorphism n(s) as a product of an appropriate 
power of a and /3. Note that e is the identity automorphism, 
and the product notation assumes the rightmost factor acts 
first. 



s = so • ■ 

7r(s) 


S 4 


40123 
a 


34012 

a 2 


23401 
a 3 


12340 

a 4 


01234 
a 5 = e 


s = so- ■ 

7r(s) 


S4 


04321 



10432 
a/3 


21043 

a 2 /3 


32104 

a 3 /3 


43210 

a 4 /3 



integer strings s = sq ■ ■ ■ S4 that result from the ten CBS 
\sb) that span the final ground-state subspace (see Sec- 
tions QTB] and [Ell . Each integer string s fixes the bottom 
row of a permutation ir(s) (see eq. (j26| ) which is a graph 
automorphism of C5 . The demonstration of this is iden- 
tical to the demonstration given for C4 and so will not be 
repeated here. Just as for Aut(Ci), the graph automor- 
phisms in Table [TTI are all the elements of Aut(C§) which 
is seen to have order 10. Note that this is the same as 
the order of the dihedral group D 5 . 

We now show that Aut(C^) is isomorphic to the dihe- 
dral group Dq by showing that the graph automorphisms 
a = tt(40123) and (3 = tt(04321) generate Aut(C 5 ), and 
satisfy the generator relations (Eq. (|27p) for D 5 . The 
second row of Table [TT] establishes that a and (3 are the 



10 



generators of Aut(C^) as it shows that each element of 
Aut(C§) is a product of an appropriate power of a and 
/3, and that all possible products of powers of a and 
(3 appear in that row. Following the discussion for C4, 
it is a simple matter to show that a corresponds to a 
72° clockwise rotation of C5. Thus 5 applications of 
a rotates C5 by 360° which leaves it invariant. Thus 
a 5 = e which is the first of the generator relations in 
Eq. (|2"71) . Similarly, (3 can be shown to correspond to a 
reflection of C5 about a vertical axis passing through 
vertex in Figure Q2] Thus two applications of (3 
leave C5 invariant. Thus j3 2 — e which is the second 
generator relation in Eq. (l27|) . Finally, using Table [Til 
direct calculation as in Eqs. ([35]) and ([36)) shows that 
a/3 = f3a 4 which establishes the final generator relation 
in Eq. (|2"T)) . We see that a and /3 generate Aut(Cs) and 
satisfy the generator relations for D5 and so generate 
a 10 clement group isomorphic to D§. In summary, 
we have shown that the GI AQA found all ten graph 
automorphisms of C5, and that the group formed from 
these automorphisms is isomorphic to the dihedral group 
D5 which is the correct automorphism group for C5. 

N = 6 : The cycle graph Cq appears in Figure [14] 
It has 6 vertices and 6 edges; degree sequence 





A = 



(38) 



FIG. 14: Cycle graph C 6 

{2, 2, 2, 2, 2, 2}; and adjacency matrix 

/0 1 1\ 
10 10 
10 10 
10 10 
10 1 

\l 1 0/ 



Numerical simulation of the GI AQA applied to the 
GI instance with G = Cq and C = G found a vanish- 
ing final ground-state energy E gs = 0. The GI AQA 
thus correctly identifies C$ as being isomorphic to it- 
self. The simulation also found that the final ground- 
state subspace is 12-fold degenerate. Table Hill lists the 



TABLE III: Automorphism group Aut{C&) of the cycle graph 
C(, as found by the GI AQA. The first and third rows list the 
integer strings s = so • • • S5 determined by the labels of the 
twelve CBS \sb) that span the final ground-state subspace 
(see text). Each string s determines a graph automorphism 
7r(s) via Eq. (|26[) . The second and fourth rows associate 
each integer string s in the first and third rows with a graph 
automorphism ir(s). They also identify the two strings that 
give rise to the graph automorphisms a and j3 that generate 
Aut(Ce), and write each graph automorphism ir(s) as a 
product of an appropriate power of a and /3. Note that 
e is the identity automorphism, and the product notation 
assumes the rightmost factor acts first. 



= «o • ■ • S5 

7r(s) 



s — so-- 

7r(s) 



*r, 



501234 450123 345012 234501 123450 012345 



105432 210543 321054 432105 543210 054321 

/3 a/3 a 2 /3 a 3 /3 a A /3 a 5 /3 



integer strings s = Sq ■ ■ ■ s 5 that result from the twelve 
CBS \sb) that span the final ground-state subspace (see 
Sections III Bl and lIII[) . Each integer string s fixes the bot- 
tom row of a permutation tt(s) (see eq. (|26])) which is a 
graph automorphism of Cq. The demonstration of this 
is identical to the demonstration given for C4 and so will 
not be repeated here. Just as for Aut{Ci), the graph au- 
tomorphisms in Table HTT1 are all the elements of Aut(Ce) 
which is seen to have order 12. Note that this is the same 
as the order of the dihedral group Dq . 

We now show that Aut(Co) is isomorphic to the dihe- 
dral group Dq by showing that the graph automorphisms 
a = tt(501234) and (3 = tt(105432) generate Aut(C 6 ), 
and satisfy the generator relations (Eq. (1271) ) for D§. 
The second and fourth rows of Table Hill establish that a 
and j3 are the generators of Aut(Ce) as they show that 
each element of Aut(Ce) is a product of an appropriate 
power of a and /?, and that all possible products of 
powers of a and j3 appear in these two rows. Following 
the discussion for C4, it is a simple matter to show 
that a corresponds to a 60° clockwise rotation of Cq. 
Thus 6 applications of a rotates Cg by 360° which 
leaves it invariant. Thus a e — e which is the first of 
the generator relations in Eq. ([2"T]) . Similarly, j3 can 
be shown to correspond to a reflection of C§ about a 
vertical axis that bisects C§. Thus two applications of 
j3 leave Cq invariant. Thus /3 2 = e which is the second 
generator relation in Eq. (|2"T)l . Finally, using Table IIII1 
direct calculation as in Eqs. (|3"5l) and (|3^1) shows that 
a/3 = /3a 5 which establishes the final generator relation 
in Eq. ([27]) . We see that a and f3 generate Aut(C§) and 
satisfy the generator relations for Dq and so generate 
a 12 element group isomorphic to Dq. In summary, we 
have shown that the GI AQA found all twelve graph 
automorphisms of Cq, and that the group formed from 
these automorphisms is isomorphic to the dihedral group 
Dq which is the correct automorphism group for Cq. 



11 



N = 7 : The cycle graph C7 appears in Figure [TS] 
It has 7 vertices and 7 edges; degree sequence 






(39) 



C 7 



FIG. 15: Cycle graph C 7 . 
{2, 2, 2, 2, 2, 2, 2}; and adjacency matrix 

/o 1 0000 l\ 

10 10 

10 10 

10 10 

10 10 

10 1 

\1 00 10/ 



Numerical simulation of the GI AQA applied to the 
GI instance with G = C7 and G' = G found a vanishing 
final ground-state energy E gs — 0. The GI AQA thus 
correctly identifies C7 as being isomorphic to itself. The 
simulation also found that the final ground-state sub- 
space is 14-fold degenerate. Table IIVI lists the integer 
strings s = Sq ■ ■ ■ sq that result from the fourteen CBS 
\sb) that span the final ground-state subspace (see Sec- 
tions QTB] and Qll|). Each integer string s fixes the bottom 
row of a permutation 7r(s) (see eq. (|26[) ) which is a graph 
automorphism of 67. The demonstration of this is iden- 
tical to the demonstration given for C4 and so will not 
be repeated here. Just as for Aut(C4), the graph auto- 
morphisms in Table llVl are all the elements of Aut(Ci) 
which is seen to have order 14. Note that this is the same 
as the order of the dihedral group Dj. 

We now show that Aut(C*i) is isomorphic to the dihe- 
dral group D-j by showing that the graph automorphisms 
a = tt(6012345) and /3 = tt(0654321) generate Aut(C 7 ), 
and satisfy the generator relations (Eq. (|27I0 for D 7 . 
The second and fourth rows of Table HVl establish that a 
and j3 are the generators of Aut(Ci) as they show that 
each element of Aut{C 7 ) is a product of an appropriate 
power of a and /3, and that all possible products of 
powers of a and (3 appear in these two rows. Following 
the discussion for C4, it is a simple matter to show that 
a corresponds to a 27r/7 radian clockwise rotation of 
C7. Thus 7 applications of a rotates C-j by 360° which 



TABLE IV: Automorphism group Aut(C-j) of the cycle graph 
CV as found by the GI AQA. The odd rows list the integer 
strings s = so • ■ ■ se determined by the labels of the fourteen 
CBS I si,) that span the final ground-state subspace (see 
text). Each string s determines a graph automorphism n(s) 
via Eq. (|26[) . Each even row associates each integer string 
s in the odd row preceding it with a graph automorphism 
7r(s). They also identify the two strings that give rise to 
the graph automorphisms a and /3 that generate Aut(C7), 
and write each graph automorphism n(s) as a product 
of an appropriate power of a and ft. Note that e is the 
identity automorphism, and the product notation assumes 
the rightmost factor acts first. 



s = s • ■ 


se 


6012345 


5601234 


4560123 


3456012 


7r(s) 




Q 


a 2 


a 3 


a 4 


s = So • • 


S6 


2345601 


1234560 


0123456 




7r(s) 




Q 5 


a 6 


7 
a — e 




s = So • • 


S6 


0654321 


1065432 


2106543 


3210654 


7r(s) 




P 


a/3 


a 2 /? 


a 3 /3 


s = So • ■ 


S6 


4321065 


5432106 


6543210 




7r(s) 




a 4 /3 


Q 5 /3 


a 6 /? 





leaves it invariant. Thus a 7 — e which is the first of 
the generator relations in Eq. (|2T|) . Similarly, f3 can 
be shown to correspond to a reflection of C-j about a 
vertical axis that passes through vertex in Figure [131 
Thus two applications of j3 leave C-j invariant. Thus 
fP = e which is the second generator relation in Eq. (|27[) . 
Finally, using Table HVl direct calculation as in Eqs. ([33)1 
and (|36|) shows that a/3 = (3a 6 which establishes the 
final generator relation in Eq. (J2"T)) . We see that a and 
(3 generate Aut(Cr) and satisfy the generator relations 
for D7 and so generate a 14 element group isomorphic 
to Dj. In summary, we have shown that the GI AQA 
found all fourteen graph automorphisms of C7, and 
that the group formed from these automorphisms is 
isomorphic to the dihedral group Dj which is the correct 
automorphism group for C7. 



2. Grid graph C?2,3 

The grid graph G2.3 appears in Figure [TBI It has 6 
vertices and 7 edges; degree sequence {3, 3, 2, 2, 2, 2}; and 
adjacency matrix 



/0 1 1 0\ 
10 10 
10 110 
110 1 
10 1 

\ooo 1 1 0/ 



.4 



(40) 



12 




G 2 



FIG. 16: Grid graph G 2 , 3 . 



G2,3 has two reflection symmetries. The first is reflection 
about the horizontal line passing through vertices 2 and 
3, and the second is reflection about a vertical line that 
bisects G2.3. Let a and (3 denote the permutation of 
vertices that these reflections produce. It follows from 
their definition that a = tt(452301) and /3 = vr(103254). 
They are the generators of the automorphism group of 
G*2,3: Aut(G2,3) = (a, (3). As reflections, they satisfy 
a 2 = f3 2 = e. Applying these reflections successively 
gives a(3 = 7r(543210) which is a fourth symmetry of 
£2,3. Applying these reflections in the reverse order, it 
is easy to verify that a (3 — (3a. Thus Aut(G2,z) is a 4 
element Abelian group generated by the reflections a and 
/3. Let us now compare this with the results found by the 
GI AQA. 

Numerical simulation of the GI AQA applied to 
the GI instance with G = G2.3 and G' — G found a 
vanishing final ground-state energy E gs = 0. The GI 
AQA thus correctly identifies G2.3 as being isomorphic 
to itself. The simulation also found that the final 
ground-state subspace is 4- fold degenerate. Table IV1 lists 



bottom row of a permutation ir(s) (see eq. ([26]) ) which 
is a graph automorphism of G2,3- The demonstration 
of this is identical to the demonstration given for G4 
and so will not be repeated here. Notice that the GI 
AQA found the graph automorphisms 7r(452301) and 
7r(103254) which implement the reflection symmetries of 
G2,3 described above. As they implement reflections, it 
follows that a 2 — (3 2 = e. The GI AQA also found the 
graph automorphism 7r(543210) which is the composite 
symmetry af3 described above. Using Table [V] direct 
calculation as in Eqs. ([33)) and (f31)|) shows that a (3 = /3a. 
Thus the GI AQA correctly found the two generators a 
and (3 of Aut{G2,2,)] correctly determined all 4 elements 
of Aut(G2.3); and correctly determined that Aut(G2,z) 
is an Abelian group. In summary, the GI AQA correctly 
determined the automorphism group of G2,3- 



3. Wheel graph W7 
The wheel graph W7 appears in Figure [17] It has 7 




W 7 



TABLE V: Automorphism group Aut{G2,'i) of the grid 
graph G2.3 as found by the GI AQA. The first row lists 
the integer strings s — so ■ ■ ■ s$ determined by the labels of 
the four CBS \sb) that span the final ground-state subspace 
(see text). Each string s determines a graph automorphism 
7r(s) via Eq. (|26[) . The second row associates each integer 
string s in the first row with a graph automorphism n(s). 
They also identify the two strings that give rise to the graph 
automorphisms a and /3 that generate Aut(G2,z), and write 
each graph automorphism n(s) as a product of an appropriate 
power of a and /3. Note that e is the identity automorphism, 
and the product notation assumes the rightmost factor acts 
first. 



s ■ 


= 8Q---8S 


452301 


103254 


012345 


543210 


7r(s) 


Q 


P 


a 2 = j3 2 = e 


a/3 



the integer strings s = s ■ ■ ■ s 5 that result from the four 
CBS \sb) that span the final ground-state subspace (see 
Sections III Bl and IIII[) . Each integer string s fixes the 



FIG. 17: Wheel graph W 7 . 

vertices and 12 edges; degree sequence {6, 3, 3, 3, 3, 3, 3}; 
and adjacency matrix 

/o 1 1 l\ 

10 10 1 

10 10 1 

1 = 10 10 1 • (41) 

10 11 

10 10 1 

\ 1 1 1 1 1 1 J 

The automorphism group Aut{Wy) is isomorphic to the 
dihedral group D§ which, as we have seen, has 12 ele- 
ments and is generated by two graph automorphisms a 
and j3 that satisfy the generator relations in Eq. (|2"T|) with 
n = 6 and respectively, rotate W7 clockwise about vertex 
6 by 60°, and reflect W7 about a vertical axis through 
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vertex 6. The graph automorphisms of Wr thus fix ver- 
tex 6. Let us now compare this with the results found by 
the GI AQA. 

Numerical simulation of the GI AQA applied to the 
GI instance with G = Wr and G' = G found a vanish- 
ing final ground-state energy E gs = 0. The GI AQA 
thus correctly identifies Wr as being isomorphic to it- 
self. The simulation also found that the final ground- 
state subspace is 12-fold degenerate. Table IVTl lists the 



TABLE VI: Automorphism group Aut(W 7 ) of the wheel 
graph Wr as found by the GI AQA. The odd rows list 
the integer strings s — so ■ • ■ se determined by the labels 
of the twelve CBS \st) that span the final ground-state 
subspace (see text). Each string s determines a graph 
automorphism n(s) via Eq. (|26[) . Each even row associates 
each integer string s in the odd row preceding it with a 
graph automorphism ir(s). They also identify the two strings 
that give rise to the graph automorphisms a and j3 that 
generate Aut(Wr), and write each graph automorphism 7r(s) 
as a product of an appropriate power of a and j3. Note that 
e is the identity automorphism, and the product notation 
assumes the rightmost factor acts first. 



s = So ■ • 


S6 


5012346 


4501236 


3450126 


7r(s) 




a 


a 2 


a 3 


s = So ■ • 


S6 


2345016 


1234506 


0123456 


7t(s) 




4 

a 


Q 5 


a 6 =e 


s = so ■ • 


S6 


1054326 


2105436 


3210546 


7t(s) 







afi 


a 2 /3 


s — so ■ ■ 


S6 


4321056 


5432106 


0543216 


7t(s) 




a 3 /3 


Q 4 /3 


a 5 /3 



integer strings s = sq ■ ■ ■ sq that result from the twelve 
CBS \si,) that span the final ground-state subspace (see 
Sections III Bl and llllj) . Each integer string s fixes the bot- 
tom row of a permutation tt(s) (see eq. (|26l) ) which is a 
graph automorphism of Wr- The demonstration of this 
is identical to the demonstration given for C4 and so will 
not be repeated here. Just as for Aut(C4), the graph au- 
tomorphisms in Table WT\ are all the elements of Aut(Wr) 
which is seen to have order 12. Note that this is the same 
as the order of the dihedral group Dq. 

We now show that Aut(Wr) is isomorphic to the dihe- 
dral group Dq by showing that the graph automorphisms 
a = tt(5012346) and (3 = tt(1054326) generate Aut(W 7 ), 
and satisfy the generator relations (Eq. (|2"T)) ) for Dq. The 
second and fourth rows of Table IVII establish that a and 
j3 are the generators of Aut(Wr) as they show that each 
element of Aut(Wr) is a product of an appropriate power 
of a and /3, and that all possible products of powers of a 
and P appear in these two rows. Following the discussion 
for C4, it is a simple matter to show that a corresponds 
to a 60° clockwise rotation about vertex 6 of Wr- Thus 
6 applications of a rotates Wr by 360° which leaves it 
invariant. Thus a 6 = e which is the first of the generator 
relations in Eq. (|27[) . Similarly, /3 can be shown to corre- 



spond to a reflection of Wr about a vertical axis passing 
through vertex 6 of Wr- Thus two applications of /3 leave 
Wr invariant. Thus /3 2 — e which is the second generator 
relation in Eq. (f2~T|) . Finally, using Table IVTl direct cal- 
culation as in Eqs. (1351) and f|36[) shows that a/3 = /3cv 5 
which establishes the final generator relation in Eq. (|27l) . 
We see that a and (3 generate Aut(Wr) and satisfy the 
generator relations for Dq and so generate a 12 element 
group isomorphic to Dq. In summary, we have shown 
that the GI AQA found all twelve graph automorphisms 
of Wr, and that the group formed from these automor- 
phisms is isomorphic to the dihedral group Dq which is 
the correct automorphism group for Wr- 



V. EXPERIMENTAL IMPLEMENTATION 

In this section we express the GI problem Hamiltonian 
Hp in a form more suitable for experimental implemen- 
tation. We saw in Section lll Cl that the eigenvalues of Hp 
are given by the cost function C(s) which is reproduced 
here for convenience: 



with 



C(s) = C 1 (s) + C 2 (s) + C 3 (s), 



7V-1 M 
N-2 N-l 

C ^ = E E s >i«> 

C 3 (s) = \\a(s)Aa T (s)-A%. 



(42) 

(43) 
(44) 

(45) 



Here s — so ■ ■ ■ sat_i is the integer string derived from the 
binary string s b = (s •• -Sy-i) • •• (s(jv-i)i/ ' • -snu-i) 
via Eqs. © and Q with U = [log 2 N] . The matrix 
a(s) was defined in Eq. ((5J) as 



<Ji,j{s) 



0, if Sj > N - 1 
&,„, if < Si < N-l. 



(46) 



Note that <Jij (s) can be written more compactly as 



u-i 



cnj(s) 



n (1 - *. 



(47) 



=JV 



We see from Eqs. (|52 ]) -pg ]) and (gTJ) that the s- 
dependence of C(s) enters through the Kronecker deltas. 
This type of s-dependence is not well-suited for experi- 
mental implementation and so our task is to find a more 
convenient form for the Kronecker delta. 

We begin with 5 a ^ in the case where a, b £ {0, 1}. Here 
we write 



Sa.b — (0 + b - iy 



(a^b) 

1 (a = b), 



(48) 
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which can be checked by inserting values for a and b. 
Now consider 5 s .k when s and fc are U-hit integers. The 
binary decompositions of s and k are 



(7-1 



i=0 

E^r. 



(49) 
(50) 



i=0 



For s and k to be equal, all their corresponding bits must 
be equal. Thus we can write 



u-x 



n *«.* 

n(^-i) 2 =(; j alu= "; 

f = o [ (some s, jt 



¥=h 



(51) 



Eq. (|5ip allows each Kronecker delta appearing in C(s) 
to be converted to a polynomial in the components of 
the integer string s. From Eqs. (|4"3"f and (14~41 we see that 
Ci(s) and C2(s) are each 2U- local. For C?3(s) we must 
write cr(s)j4(T T (s) in a form that makes the Kronecker 
deltas explicit. Using Eq. (|47|) . we have 

JV-l 
<7(s)ykr T (s) = E cr/ i (s)AyCr m: ,-(s) 

iV-1 r £/-l ^ 

= e ^ n (i-^.^u^ 

i,j=0 L a=N ) 

U-X 

6m, Sj J] C 1 " **./»)>• ( 52 ) 



/3=N 



Inserting Eq. (j5"Tj) into Eq. (15^1) . we see that a(s)Aa T (s) — 
A' is 4U(U - N + l)-local. If we use the Li-norm (L 2 - 
norm) in Eq. ((45]), then we see that C 3 (s) is 4U(U — N + 
l)-local (8C/(C7 - N + l)-local). To implement C(s) on 
the D-Wave quantum annealing device, C(s) must be re- 
duced to 2-local, and the resulting qubit-qubit couplings 
must match the hardware's Chimera coupling-graph. A 
procedure for carrying out this reduction is detailed in 
Ref. pjj. 



VI. SUBGRAPH ISOMORPHISM PROBLEM 

An instance of the SubGraph Isomorphism (SGI) prob- 
lem consists of an TV-vertex graph G and an n-vertex 
graph H with n < N. The question to be answered 
is whether G contains an n-vertex subgraph that is iso- 
morphic to H . The SGI problem is known to be NP- 
Complete [2(| and so is believe to be more difficult to 
solve than the GI problem. Here we show how an instance 
of SGI can be converted into an instance of a COP whose 



cost function has a minimum value that vanishes when 
G contains a subgraph isomorphic to H, and is greater 
than zero otherwise. The SGI cost function will be seen 
to be a natural generalization of the GI cost function 
given in Eqs. ([9])— f[T2]>. The SGI COP can then be solved 
using adiabatic quantum evolution as was done for the 
GI problem. 

As with the GI problem, we would like to determine 
whether there exists an isomorphism 7r of G that pro- 
duces a new graph ir(G) that contains H as a subgraph. 
Just as with the GI problem, we consider linear maps 
a(s) (see Eqs. ©-©) that transform the adjacency ma- 
trix A of G to A(s) = a{s)Aa T (s). We then search n(G) 
to determine whether there is a subset of n vertices that 
yields a subgraph that is equal to H . 

To begin the process of converting an SGI instance into 
an instance of a COP, let: (i) a label all the ( ) ways 
of choosing n vertices from the N vertices in n(G); and 
(ii) \i) (\o.i}) be an n-component (iV-component) vector 
whose i-th (a^-th) component is 1, and all other compo- 
nents are 0. Thus a labels the choice (aoi ■ • ■ > O-n— l) of n 
vertices out of the N vertices of 7r(G). We now show that 
an n x N matrix P a can be used to form annxn matrix 
A a (s) whose matrix elements are the matrix elements of 
A(s) associated with the n vertices appearing in a. To 
that purpose, define 



p q = Enx^i. 



i=0 



A a (s) = P a A{s)Pl , 



where 



N-X 



A(s)= E Am(s)\l){r 



(53) 



(54) 



(55) 



l,m=0 

It follows from these definitions that 

n-X N-X n—X 

a q (s) = EN)(^I E 4m(«)IOHX>'-WI 

i.=0 Lm=0 j—0 

n-X 

= E^oowo'i- (56) 

ij'=0 

Thus the matrix elements (A a )ij (s) are precisely the ma- 
trix elements A ai , a .(s) associated with all possible pairs 
of vertices drawn from (a.o, . . . , a n _i). The matrix A a (s) 
is thus the adjacency matrix for the subgraph g a com- 
posed of the vertices appearing in a, along with all the 
edges in 7r(G) that join them. With A a (s), we can, as in 
the GI problem, test whether g a is equal to H by check- 
ing whether ||.4 a (s) — A'\\i vanishes or not. Here A' is 
the adjacency matrix of the graph H, and ||C||i is the 
Li-norm of O. 

We can now define a cost function whose minimum 
value vanishes if and only if G contains a subgraph g 
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that is isomorphic to H. Because the transformation 
cr(s) must be a permutation matrix when g is isomorphic 
to H, we again introduce the penalty functions Ci(s) 
and Cz(s) used in the GI COP to penalize those integer 
strings s that produce a a(s) that is not a permutation 
matrix: 



N-l M 

C ^*) = E E *«*.<• 

i=0 a=JV 
JV-2 JV-1 

^w = £ E *«»* 

i=0 j=i+l 



(57) 
(58) 



The final penalty function Cs(s) for the SGI problem 
generalizes the one used in the GI problem. It is defined 
to be 



(?) 

C 3 (s) = l[\\A a (s)-A'\ 



(59) 



Q = l 



/N\ 



where the product is over all ( j ways of choosing n 
vertices out of N vertices. Note that C 3 (s) vanishes if and 
only if G contains an n- vertex subgraph isomorphic to H. 
This follows since, if G contains an n-vertex subgraph g 
isomorphic to H , there exists a permutation n(s) of G 
that has an n-vertex subgraph that is equal to H. Thus, 
for this s, there is a choice a of n vertices that gives 
a subgraph for which A a (s) — A' = 0. It follows from 
Eq. flS§) that C 3 (s) = 0. On the other hand, if C 3 (s) = 0, 
it follows that at least one of the factor on the RHS of 
Eq. (|5^| vanishes. Thus there is a choice a of n vertices 
for which ||A*(s) - A'\\i = 0. Thus A a (s) = A, and so 
G has a subgraph isomorphic to H. 

The cost function for the SGI problem is now defined 
to be 



(60) 



C(s)=C l (s) + C 2 (s) + C 3 (s), 



where C\(s), C2(s), and C 3 (s) are defined in Eqs. (|5"T)) - 
(15^1) . This gives rise to the following COP: 

SubGraph Isomorphism COP: Given an TV-vertex 
graph G and an n- vertex graph H with n < N, 
and the associated cost function C(s) defined in 
Eq. (|60|) . find an integer string s* that minimizes 
C(s). 

By construction: (i) C(s*) = if and only if G contains a 
subgraph isomorphic to H, and er(s*) is the permutation 
matrix that transforms G into a graph tt(G) that has H 
as a subgraph; and (ii) C(sA > otherwise. 



As with the GI COP, the SGI COP can be solved us- 
ing adiabatic quantum evolution. The quantum algo- 
rithm for SGI begins by preparing the L — N\\og 2 N~\ 



qubit register in the ground-state of Hi (see Eq. (fTSj) ) and 
then driving the qubit register dynamics using the time- 
dependent Hamiltonian H{t) = (1 - t/T)H l + (t/T)H P . 
Here the problem Hamiltonian Hp is defined to be diago- 
nal in the computational basis \sb) and to have associated 
eigenvalues C(s), where s is found from Sb according to 
Eqs. ([5]) -0. At the end of the evolution the qubits are 
measured in the computational basis. The outcome is 
the bit string s£ so that the final state of the register 
is \sl) and its energy is C(s*), where s* is the integer 
string derived from sj*. In the adiabatic limit, C(s*) will 
be the ground-state energy, and if C(s*) = ( > 0) the 
algorithm concludes that G contains (does not contain) 
a subgraph isomorphic to H . In the case where G does 
contain a subgraph isomorphic to H, the algorithm also 
returns the permutation 7r» = 7r(s*) that converts G to 
the graph 7r*(G) that contains H as a subgraph. Note 
that any real application of AQE will only be approxi- 
mately adiabatic. Thus the probability that the final en- 
ergy C{s*) will be the ground-state energy will be 1 — e. 
In this case the SGI quantum algorithm must be run 
k ~ C(ln(l — 6) fine) times so that, with probability 
5 > 1 — e, at least one of the measurements will return 
the ground-state energy. We can make S arbitrarily close 
to 1 by choosing k sufficiently large. 



VII. SUMMARY 

In this paper we have presented a quantum algo- 
rithm that solves arbitrary instances of the Graph Iso- 
morphism problem. We numerically simulated the algo- 
rithm's quantum dynamics and showed that it correctly: 
(i) distinguished non-isomorphic graphs; (ii) recognized 
isomorphic graphs; and (iii) determined the automor- 
phism group of a given graph. We also discussed the 
quantum algorithm's experimental implementation, and 
showed how it can be generalized to give a quantum algo- 
rithm that solves arbitrary instances of the NP-Complete 
SubGraph Isomorphism problem. 
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